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ANALYSIS  OF  THE  PUFFS  TFIACKER  WOX-IA 


AS  A  GENERAL  NONLINEAR  FILTER 


Prepared  by: 

C.  Nicholas  Pryor 


ABSTRACT:  A  statistical  analysis  is  made  of  the  PUFFS  Target 

Tracker  WOX-IA.  This  device  uses  a  nonlinear  tracking  scheme 
such  that  the  estimated  target  position  is  changed  in  discrete 
steps  of  fixed  magnitude,  dependent  only  upon  the  sign  of 
observed  error  at  regular  sampling  intervals. 

The  theoretical  performance  characteristics  of  the  tracker 
are  determined  from  the  statistical  analysis.  Equations  are 
developed  for  the  output  standard  deviation  of  tracker  position, 
the  lag  in  tracker  position  due  to  bearing  rates,  and  the 
equivalent  integration  time  of  the  tracker  as  functions  of  the 
input  standard  deviation  and  the  tracker  step  size.  These 
equations  are  expressed  in  a  sufficiently  general  form  so  that 
they  can  be  readily  applied  to  other  systems.  Refinements  are 
also  discussed  which  eliminate  the  bearing  lag  problem,  and 
stability  of  the  rate-correction  system  is  investigated. 

Finally  experimental  verification  is  presented  for  some  of 
the  results,  and  the  qualitative  operation  of  the  tracker  is 
discussed  as  an  adaptive  system  interchanging  bandwidth  and 
output  noise. 

Other  applications  for  nonlinear  systems  of  this  type  are 
mentioned,  including  tracking  and  post-filtering  problems  in 
sonar  or  radar,  and  filtering  problems  in  navigation  systems. 
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This  report  is  a  statistical  and  information  theory  analysis 
of  the  performance  of  an  automatic  target  tracker  developed  at 
the  Naval  Ordnance  Laboratory  in  which  it  is  shown  that  the 
tracker  is  a  special  case  of  a  group  of  generalized  nonlinear 
filters.  This  work  was  funded  under  Task  No.  TfF008-04-004-RUDC- 
3C-150,  PUFFS  Technical  Direction.  This  report  will  be  of  special 
interest  to  those  concerned  with  the  sonar  automatic  tracking 
problem  and  to  specialists  in  signal  processing.  However,  the 
theory  should  also  be  useful  in  radar  and  navigation  systems  and 
has  possible  applications  in  the  fields  of  communications,  IFF 
and  detection. 


W.  D.  COLEMAN 
Captain,  USN 
Commander 
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ANALYSIS  OF  THE  PUFFS  TRACKER  WOX-IA 
AS  A  GENERAL  NONLINEAR  FILTER 

INTRODUCTION 

The  automatic  tracker^  developed  for  the  PUFFS  sonar  system 
is  used  to  follow  the  position  of  the  correlogram  in  the  PUFFS 
display,  and  to  smooth  out  the  major  fluctuations  of  the  correl¬ 
ogram  which  occur  due  to  noise  in  the  environment  and  finite 
integration  time.  The  operation  of  the  tracker  is  basically  as 
follows:  Each  time  a  correlogram  is.  generated  (once  every 

25  milliseconds)  the  tracker  makes  an  area  measurement  to  deter¬ 
mine  vv^ether  the  center  of  the  correlogram  is  to  the  right  or  to 
the  left  of  the  last  tracker  estimate  of  the  position.  On  the 
basis  of  this  measurement,  the  tracker  then  changes  its  estimate 
by  a  fixed  amount  in  the  direction  toward  the  observed  center  of 
the  correlogram.  This  is  the  nonlinear  feature  of  the  tracker; 
a  linear  system  would  move  a  distance  proportional  to  the 
observed  error,  vdiereas  the  tracker  makes  a  correction  of 
constant  magnitude  regardless  of  the  magnitude  of  the  error. 

A  post  integrator  for  the  PUFFS  system  has  also  been  proposed 

in  which  samples  of  the  observed 'correlogram  are  compared  with 

...  stored -estimates  of  the  correlogram  at  a  number  of  corresponding 

delay  times.*  Then  the  estimates  are  increased  or- decreased  a 

fixed  amount  depending  on* whether  the  observed  sample  is  larger 

or  smaller  than  the  stored  values  for  the  given  delay.  It  may  be 

seen  that  both  these  devices  are  similar  in  principle  and  may.  be 

generalized  into  a  class  of  nonlinear  filters.  Such  a  filter  • 

would  repeatedly  compare  its  stored  estimate  of  the  input 

variable  with  samples  of  the  input  and  then  change  the  stored 

value  at  each  comparison'  by  ci  fixed  step  size  in  the  direction 

of  the  observed  input  value.  The  analysis  will  therefore  be  in 

terms  of  this  generalized  filter,  and  the  results  will  apply 

equally  well  to  the  tracker,  pos-t  integrator,  or  any  other  device 

operating  in  a  similar  manner.  Although  the  analysis  here  is 

based  on  a  sampled  system  with  quantized  output  values,  the 

results  may  be  extended  to  continuous  systems  with  continuous' 

output  distributions.  The  sampled  approach  was  chosen  because 

it  affords  more  direct  analysis  without  discussion  *of  frequency 

response  of  the  .elements,  and  because  it  corresponds  to  the 

cases  of  interest  in  the  PUFFS  system. 

■ 


1 - ^ - 

Munson,  Cope,  Sabelhaus,  "Target  Tracker  ¥K)X-1A”  ,  NAVORD 

Report  6799  (in  preparation). 
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THEORY  OF  OPERATION 

One  j{  the  first  items  of  interest  in  discussing  the  opera¬ 
tion  of  such  a  filter  is  the  standard  deviation  of  the  noise 
output  of  the  system,  when  white  noise  of  a  known  standard 
deviation  is  applied  to  the  input.  We  may  determine  this  by 
computing  the  probability  of  finding  the  output  of  the  filter  in 
each  of  its  possible  output  positions  or  states.  If  the  filter 
output  is  in  its  kth  state  so  that  the  output  value  Xq  of  the 
filter'is  Xq  =  Ak,  where  A  is  the  size  of  the  steps  and  the  0th 
state  is  defined  to  be  at  x©  =  0,  then  the ^probability  of  the 
input  sample  being'  less  than  this  value  is^Jby  definition  F(xo) 
or  F{Ak)  where  F(x)  is  the  cumulative  distribution  function  of 
Xq.  This^'is  the  probability  that  the  filter  will  decrease  its 
output  value  by  the  step  A.  Similarly  the  filter  will  increase 
its  output  with  a  probability  l-F(Ak)  or  the  probability  that 
the  input  sample  is  larger  than  Ak.  If  we  assume  that  the 
system  is  in  its  steady -.state ,  then  the  probability  of  moving 
from  its  kth  state  to  its  k  +  Ist  state  must  equal  the  proba¬ 
bility  of  moving  from  the  k  +  1st  state  into  its  kfh  state.  The 
probability .  of  moving  from  k  to  k  +'  1  -is  equal  to  the  proba¬ 
bility  of  being  in  .the  kth  state  P|(  times  the  probability  of 
increasing  its  value  or- 1  -  F(Ak) .  The  probability 'of  moving 
from  k  +  1  to  k  is  Pj^  1  times  the  probability  of  decreasing 
its  value  or  F(A(k  +  1)).  Equating  these  we\have : 

Pk  (1  -  F(Ak))  =  Pk  +  1  F(A(k  +  D)  or 

This  is  the  basic  equation  relating  the  state  probabilities  for 
the  filter.  By  iterating  equation  (1)  we  can  find  the  proba¬ 
bility  of  any  state  k  in  terms* of  Pq  since 

P  =  P  1  7,f(°)  .  . 1  -  -  1)) 

o  FIA)  F(2A)  t(Ak) 


This  can  be  written  in  the  form: 


TT'  F(Aj) 

j  =  l 


If  we  assume  the  input  to  be  white  noise  having  a  Gaussian  dis¬ 
tribution  with  a  mean  of  zero  and  a  standard  deviation  crj  then 
for  values  of  x  small  compa’red  to  crj^  the  distribution  function 
F.(x)  becomes  approximately 
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Inserting  this  in  equation  (2)  we  obtain 

-  yl  n) 


TT 


p,  =  p 


JT-O- 


=  p 


J=o 


(l  -  aj) 


^  A  . y||i)  °  ^ 

j  =  l'  tt  ai/  ^  ^ 


(3) 


where 


a  -  A^ 
✓  TT  Oi 


Since 


TT  (l+  a j)  =  V  +  a j)=  1  +  a  j  +  {  )  +=  1  +  a  +  (  ) 

j=o  j  =  i''  •  j=l 


we  can  evaluate  equation  (3)  as: 

=  P 

o 


P,  ^  P 


and  substituting  for  a  we  obtain  the  relative  state  probabilities 
or: 


P  =  P  (1  -  ~  k^) 

^.k  ^^O  ^  ■/  TT  csi  ^  ' 


(4) 


Since  the  state  P^  has  an  output  value  Ak,  we  can  write  the  pro¬ 
babilities  of  the  output  values  as 


P(Ak)  =  P(o)  (1  - 


'2  A  ^2) 


V  Oi 


Finally  letting  Xq  equal  Ak  we  have  the . probabi lity  of  the  out¬ 
put  value  Xq  where  Xq  is  restricted  to  integfal  multiples  of  A. 


P(xo)  =.P(o)  (1  -/igj 


(5) 


Although  equation  (5)  has  a  parabolic  form,  it  has  been  obtained 
from  truncation  of  a  Taylor  series,  and  it  would  be  reasonable 
to  suspect  that  the  true  form  of  P(xo)  is  approximately  Gaussian. 
This  is  also  borne  out  by  experiment  and  by  computation  of  P(xo) 
directly  from  equation  (1).  Using  this  assumption,  we  may 
compare  equation  (5)  with  the  Taylor  expansion  of  the  Gaussian 
distribution : 


1 
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2a. 


o  = 


y2iT 


(1  - 


2a_2 


Equating  the  coefficients  of  the  squared  terms  in  the  series 
expansion  we  obtain: 

21  .2 

TT  Acri 


or 


I 


Thus  the  output  sltandard  deviation  a©  of  the  filter  has  been 
determined  in  terms  of  the  input  standard  deviation  cr^  and 
the  step  size  Ahas: 


a 

o 


*/a~[ 


(6) 


The' result  is  valid  for  cases  where  the  step  size  A  is  small 
compared  to  cfj^. 


Notice  that  the  output  standard  deviation  is  proportional  to 
the  square  root  of  the  input -standard  deviation.*  This  is  the 
striking  feature. o‘f  this  filter,  since  linear  f  liters^  would  have 
a  direct  proportionality  between  input  and  output  st^ruiard 
•  deviations."  The  significance  of  this  effect  will  be  discussed 
lat^r  in  this  report.^ 


Aside  from  the  noise  output  of  the  filter,  the  o.utput  also 
differs  from  the  ideal  input  in  that  It  lags  behind  the  input 
whenever  the  input  iias  some  rate  of  change,  p  units  per  second. 
This  occurs  because*  the  output  must  be  offset  from  the  center 
of  the  input  dens.ity  function  in  order  for  the  probability  of 
.  moving  in  the  direction  of  input  motion  to  be  signif  ica'nt  ly ' 

greater  than  the  probability  of  moving  in  the  opposite*  direction.* 
In  particular,  if  the  time  interval  bet’ween  consecutive  decisions 
is  T  seconds  and  the  rate”' is  p,  the  probability  P  of  moving  in 
the**directioTr."of  "input  motion^must  be 

=  5 ^ Is T  ■  '  (7) 

A  — wV-;—  *•*•••  • 

in  order  that  the  probability  of  moving  in  the  proper  direction 
will  exceed  the  probability  of  moving  in  the  opposite  direction 

by.^  T* and  the  net-rate  of  output  motion  will  match  the  rate  of 

input  motion  p.  The  "lag  must  be  such  that  the  areas  left  and 
right  of  the  steady  state  point  on  the  input  noise  density 
function  produce  the  above  probability.  For  small  lags  t  and  a 
standard  deviation  of  the  normally  distributed  input  signal  a^, 
the  probability  of  moving  in  the  proper  direction  (opposite  the 
lag)  is  • 
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P  =  i 

2  CTi 

Equating  equations  (7)  and  (8)  we  have 

P“  —  =  —  ^T  =  - - - 

2  2  A  *  ^  a, 


(8) 


Hence  the  lag  is  related  to  the  input  rate  and  standard  devia¬ 
tion  and  the  step)  size  by 


This  represents  a  constant  bias  error  due  to  constant  rate  of 
changes  p  in  the  input  signal. 


FILTER  WITH  RATE  CORRECTION 

• 

The  problem  of  steady  state  .lag  may  be  eliminated  by  having 
the  filter  automatically  correct  itself  for. the  observed  rate 
of  change  of  either  input  or  output.  This  may  be  done  by 
measuring  the  relative  probabilities  of  increase  and  dec.rease--" 
decisions • required  to  follow  the  input ’ signal ,  and  using  tnis 
measurement  to  make  adcTitional  corrections  ‘in  the  output  rate  • 
to  compensate  for  the  lag.  The . particular  solution  chosen  in 
the  tracker2  was  to  store  a  rate  estimate  as  well  as  a  position 
estimate  in  the  tracker  and  use  the  same  error  signals  to  change 
both  estimates;  Thus  an  observed  error  in  the  negative  direc-  , 
tion  will  increase  the  position  e-stimat-e -by-an- amount  A  and’ 
increase  the  rate  estimate  by  an  amount '9.  This  rate  estimate 
is  inserted  in  the  position  estimate  along  with  the  original 
position  decisions.  Thus  in  the  steady' state  with  the  input 
changing  at  a  constant  rate,  the  rate • estimate  must  stay  con¬ 
stant  and  there  mus.t  be  as  many  increase  as  decrease  decisions 
going  into  both  rate"  and  position'  estimates .  This  means  that 

P  must  be  -x,  and  by  equation  (8)  the  steady  state  lag  must  be 
zero. 


The  filter  with  the  rate  c'orrection  added  may  be  analyzed  by 
drawing  a  logical  diagram  of  the  filter  as  shown  in  Figure  1  and 
converting  it' to  an  equivalent  linear  system.  The  on ly _ nohlineay 
part  of ’  the  system  is  the  clipper  and  for  f requencies' low  cofn- 
.  pared  to  the  sampling  frequency  we  can  approximate  it  with  a 
linear  component  as  follows; 


2 

(a  ra te -compensated  tracker  is  presently  being  developed.) 
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It- 


Assume  we  have  a  clipper  with  an  input  quasi-static  signal  Xs 
and  in  input  noise  xn  with  a  standard  devia^i-on  cr^ .  Also  assume 
the  output  of  the  -clipper  is  plus  one  if  th*e  input  is  positive 
and  minus  one  if  the  input  is  negative.  If  P(-)  is  the  proba¬ 
bility  that  the  input  +  Xg  is  negative  (and  thus  that  the 
output  is  negative),  then  the  average  output  will  be  given  by 
1,-  2P(-K  But  if  the  distribution  function  of  the  noise  is 
F(x),  then  P(-)  is  the  probability  that  Xp  ^  -  Xg  or 
p(-)  =  F(-Xs).  Thus  the  average  output-  becomes  1  -  2F(-xs). 

For  a  Gaussian  noise  input  of  standard  deviation  a  , 


F(x) 


and  the  average  output  becomes  x^ 


Thus  for  low  frequencies  and  for  signals  small  compared  to  the 
noise,  the  clipper  may  be  replaced  by  an  amplifier  whose  gain’ 
is  given  by  ^  _  /2  _1_.  For  frequencies  low  compared  to  the 

sampling  frequency  the  summers  may  also  be  replaced  by  inte¬ 
grators  with  gains  appropriate  to  the  sampling  frequency  and 
the  samplers  removed.  The  resulting  system  is  shown  in  ’Figure  2 
and  is  linear  and  continuous  in  both  amplitude  and  t-ime. 


The  response  characteristics  of  the  filter  may  now’ be 
determined  in  the  Laplace  transform  domain  simply  by  repla.cing 
the  integrators  by  1/s.  The  forward  gain  of  the  system  becomes 


TS 


H(s)  = 


[•)  and  the  closed  loop  transfer  function  is: 


—  (A  -f  -) 
TS  ^  s^ 


1  .  ^  (A  . 


GAs  +  C6 
Ts^  +  GAs  +  GG 


(10) 


It  may  be  seen  that  the  frequency  response  of  the  filter  is 
that  of  a  damped  second  order  system  where  the  resonant 
frequency  is  determined  by  the  rate  step  size  0  and  the  damping 
is  due  to  the  position  steps  A.  However  the  gain  G,  which  is  a 
function  of  the  input  standard  deviation,  is  a  factor  in  both 
resonant  frequency  and  damping  ratio.  Thus  the  parameters  A 
and  0  must  be  chosen  to  satisfy  the  desired  stability  and 
frequency  response  requirements  for  all  anticipated  values  of 
input  standard  deviation.  In  particular  for  critical  damping 
we  must  require  that: 


^  v/8^F  Oi  T  (11) 

at  the  largest  anticipate d  value  of  Oi ,  where  G  has  been 
replaced  by  /2  1  . 

y  TT  a  I 


i 
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COMPARISON  OF  LINEAR  AND  NONLINEAR  FILTER  PERFORMANCE 

t 

An  index  which  may  be  used  to  express  the  performance  of  a 
nonlinear  filter  is  the  equivalent  integration  time  of  the 
filter  compared  to  a  simple  linear  filter.  This  integration 
time  is  defined  as  the  time  constant  of  a  simple  RC  low-pass 
filter  which  .would  qive  .equivalent  performance. 

\  ' 

I  ■  ,  I 

The  most 'direct  apprjbach  for*^  finding  the  equivalent  inte¬ 
gration  time  comes  from  fthe  transfer  function  derived  for  the 
filter  in  equation  ( 10)  For  the  case  of  the  filter  without 
rate  correction  (that  is*  9  =  0),  the  transfer  function  of  the 
filter  would  be 


H(s)  =  - ■  .  . 

}  i  . “  — ■' 

•**  • 

We  may  compare  this  with  the  transfer  function  of  a  simple  RC 
integrator  with  time  constant  r,  which  is: 


H(s) 


1 

-rs  +  -1 


Comparing  like,  terms  in  the  tvr<5‘  equations ,  -and*,  remembering  that 

/2  1  •  ■ 

G  =  y  —  we  see  that  the  small-signal  transfer  function  of 

the  nonlinear  filter  is  equivalent  to  that  of  an  RC  integrator 
with  a  time  constant  given  by: 


T 


2  a 


T 


(12) 


An  independent  check  on  this-  result  may.  be  obtained  from 
equation  (9)  for  the  lag  due  to  a  constant  rate  of  change  of 
the  input.  It  may  be  shown  that  a'  simple  RC  integrator  has  a 
lag  given  by  px .  Substituting  this  for  the  lag  shown  in  equa¬ 
tion  (9)  and  solving  for  the  x  which  will  give  equiyalent-  lags, 
we  obtain  the  same  result  given  in  equation  (12). 


It  is  possible  at  this  point  to  compare  the’  noise  output 
given  in’  equation  (6)  with  the  noise  output  from  the  RC  inte¬ 
grator  with  equivalent  .time  constant  t.  For  the  RC  filter  the 
output  at  a  given  time  t  is  given,  by: 

T  ...» 


X^(t)  = 


(1  - 


z  - 

n=o 


"*e 


nT 

X 


(t  .--nT) 


(13) 


\ 
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If  we  assume  that  each  input  sample  has  a  variance 
output-  variance  may  be  computed  as  the  sum  of  the  weighted 
input  variances,  assuming  the  input  samples  are  independent 
(white  noise).  Thus: 

X  ou  X 

=  (1  -  e  ■•')2  '2-  U 
n=o 

-to  2 

=  (1  -  e  •')2  (T  2  2  e  ^ 

n=o 


=  (1  -  e 


(14) 


1  -  e 


For  sufficiently  large  y,  this  reduces  to: 

•,.2 


a  2  =  .  2  V 
O  i  oX 


=  CT 


‘^O  =  ^i  /  ^ 


2  X 

i  2'C 


( 15) 


(16) 


Substituting  t  from  equation  (12)  into  equation  (16)  we  obtain 
the  output  standard  deviation  from  the  equivalent  linear  filter 


a  = 
o 


*/t: 


(17) 


Comparing  this  to  equation  (6),  we  find  that  the  noise  output  of 
the  nonlinear  filter  is  greater  than  that  of  the  linear  filter 
with  equivalent  transfer  function  by 

*  ✓  5 

Another  method  of  expressing  the  performance  of  the  nonlinear 
f ilter  jcompared  t-o  that  of  a  linear  filter  is  to  equate  the 
output  standard  deviations  in  equations  (6)  and  (16),  and  solve 
for  the  time  constant  of  the  linear  filter  thajt  would  have 
equivalent  output  noise.  The  resultant  time  constant  is; 

-  =  /I  ^  T  •  '  (18) 


(18) 


which  is  less  than  equation  (12)  by  — .  That  is,  for  equivalent 

output  noise  the  nonlinear  filter  has  a  larger  effective  time 
constant,  and  thus  less  bandwidth,  than  the  linear  filter. 
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i 

It  is  -interesting  to  note  in  the  case  where  the  output  standard 
deviationsof  the  two  filters  are  equated  that  the  step  size  A 
used  in  the  nonlinear  filter  is  exactly  equal  to  the  expected 
magnitude  of- the  steps  in  the  linear  system,  which  is  given  by 
/2  T 

y—  —  a^.  Throwing  away  the  magnitude  and  using  only  the 

polarity  information  in' the  error  signal  resul'ts  in  a  decrease 
•in  bandwidth,  or  if  the  ?tep  size  is  adjusted  for  equal  band¬ 
width  the  output  noise  i^  increased. 

/  . 

EXPERI'MErNTTAL  VERSIFICATION  . 

I 

I-n  order”  to  provide^  an  experimental  check  on  the  analytical 
results  for  the  filter,  a  program  was  written  for  the  IBM  704 
computer  which  simulated  the  tracker  operation  and  measured  the 
“outpnt~s-bandard  deviatioh  as  a  function  of  the  step  size  A.  A 
constant ‘ input  standard  deviation  of  unity  was  used  in  the 
experiment.  The  results,  of  the  experiment,  covering  step  sizes 
from  0.1  to  0.001,  are  shown  in  Figure  3  along  with  the  theo¬ 
retical  result.  The  experimental  standard  deviation  was 
consistently  llj^  to  15%  above  the  expected  value.  This  may  be 
due  to  non-independence  of  the  input  samples  (the  input  was  a 
list  of  10,000  numbers  generated  by  a  random  number  routine.)  or 
may  be  a  result  of  the  assumptions  made  in  deriving  equation  (6). 
It  is  clear  from  Figure  3,  however,  that  the  form  of  equation  (6) 
is  correct  and  differs  at  most  by  a  small  multiplicative  constant 
from  the  correct  result.  Experiments  run  on  the  PUFFS  WOX-IA 
Tracker  have  verified  equations  (6)  and  (9)  within  the  accuracy 
to  which  the  input  standard ^ deviation  is  known  (in  fact 
equation  (9)  is' used  to  measure  the  standard  deviation  of  the 
input  to  the  tracker),  and  equations  (10)  and  (11)  have  been 
verified  through  experiment-s  on  th,e  rate  correction  circuits  for 
the""tra'cker . 


DISCISSION  AND  CONCLUSIONS 

The  analytical  work  done"  on  this  nonlinear  filter  shows  that 
the  main  feature  distinguishing  it  from  a  linear  filter  is  the 
variation  of  its  effective  integration  time  with  input  standard, 
deviation.  Thus  as  the  input  noise  increases  the  effective 
integration  time  is  increased,  with  the  result  that  the  output 
noise  only  increases  as  the  square  root  of  the  input  noise. 
System  frequency  response  is  thus  maintained  as'  long  as  the 
input  signal'  is  good  enough  to  permit  it,  and- as  the  input  noise 
increases  bandwidth  is  sacrificjed  to  maintain  a  reasonable 
output  noise  level.  The  most  direc-t  effect  of  the  loss  of  band¬ 
width  is  increased  lag,  but  as  we  have  seen  it  is  possible  to 
provide  rate*  correction  to'  eliminate  this  problem.  The  filter 
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may  then  be  considered  as  an  adaptive  system  capable  of  adjust¬ 
ing  its  bandwidth  automatically  to  reach  some  sort  of  happy 
compromise  between  bandwidth  and  output  noise  as  the  input  noise 
changes.  The  price  paid  for  this  adaptive  feature  is  a  slight 
increase  (about  23?^)  in  the  output  noise  of  the  system  if  the 
dynamic  response  is  kept  the  same,  or  a'  reduction  of  the  band¬ 
width  of  the  system  by  a  factor  of  0.66  if  output  noise  is 
maintained  at  the  same  level. 

•  • 

This  sort  of  filter  is  particularly  adaptable  to  the  PUFFS 
tracker  and  to  pther  tracking  or  post  filtering  problems  in 
sonar  or-.radar  systems,  and  to  filtering  problems  in  navigation 
systems.  Typically  in'"  t'hese  -problem,? ,  ._as  the  range  to  the  tar¬ 
get  of  the  navigational  facility  incrreaseV’t'he  noise  in  the  raw 
data  increases,  but  at!  the  same  time  the  rates  of  change  of 
input  data  decrease  ar^  the  required  bandwidth  is  t.hereby 
reduced.  Thus  a  filter  such  as  the  one  desc'ribed  malces  use  of 
the  lower  bandwidth  requirement  to  provide  additional  smoothing 
of  the  raw-data  and  considerably  extend  the /operating  range  of 
the  -entire  system.  Use  of  this  concept  in  |he  WOX-IA  PUFFS 
Tracker  has  proved  the  usefulness  of  such  filters  by  considerably 
increasing  the  tracking  capabilities  over  previous  linear  systems. 


/ 

/  . 
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LIST  OF  SYMBOLS 


Symbol 

Use 

Dimensions 

^1 

Input  Variable 

. units  ( of  input 
variable ) 

X 

0 

Output  variable 

units 

X 

s 

Signal  Component  of  Input 

units 

s 

X 

n 

Noise  Component  of  Input 

units 

^  . 

Input  Standard  Deviation  . 

units 

^0 

Output  Standard  Deviati'on 

. units 

A 

Position  Step  Size 

units 

e 

Rate  Step  Size 

units  /second 

k 

Output  State  Number 

’ dimensionless 

F(x) 

Cumulative  Distribution 
Function 

dimensionless 

P|^  or  P(x) 

Probability  of  Output  State 
or  Value 

dimensionless 

P 

Rate  of  Change  of  Input 

uni ts /second 
; 

1.  ' 

Lag  of  Output  behind  Input 

uni  ts 

G 

Gain  of  Clipper  Analogue 

—  units'^ 

V  TT 

s 

Laplace  Operator 

seconds  ^ 

H(s) 

Transfer  Function  of. ..Filter 

dimensionless 

T 

Effective  Integration  Time 

(Exponential  Time  Constant; 

seconds 

1 

T 

Sampling  (Decision)  Period 

sec  onds 
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FIG.3  OUTPUT  STANDARD  DEVIATION  VS.  STEP  SIZE 


NOL-IBM  704  SIMULATION  25  NOV.  I960 
INPUT  STANDARD  DEVIATION  =  1.00 
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